function upsilon = fn_upsilon(phi,k0)

upsilon = (phi./2).*((k0).^2);

end